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We use cosmography to present constraints on the kinematics of the Universe, without postulating 
any underlying theoretical model. To this end, we use a Monte Carlo Markov Chain analysis to 
perform comparisons to the supernova la Union 2 compilation, combined with the Hubble Space 
Telescope measurements of the Hubble constant, and the Hubble parameter datasets. We introduce 
a sixth order cosmographic parameter and show that it does not enlarge considerably the posterior 
distribution when comparing to the fifth order results. We also propose a way to construct viable 
parameter variables to be used as alternatives of the redshift z. These can overcome both the 
problems of divergence and lack of accuracy associated with the use of z. Moreover, we show that 
it is possible to improve the numerical fits by re-parameterizing the cosmological distances. In 
addition, we constrain the equation of state of the Universe as a whole by the use of cosmography. 
Thus, we derive expressions which can be directly used to fit the equation of state and the pressure 
derivatives up to fourth order. To this end, it is necessary to depart from a pure cosmographic 
analysis and to assume the Friedmann equations as valid. All our results are consistent with the 
ACDM model, although alternative fluid models, withnearly constant pressure and no cosmological 
constant, match the results accurately as well. 
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I. INTRODUCTION 

Ever since the pioneering works of two separate groups 
in 1998 cosmological observations indicate a late 
time accelerated Universe. More recently, additional evi- 
dence coming from other experiments 0, H| confirms that 
the Universe is going through an accelerated expansion. 
Thus, the existence of the acceleration is assumed to be 
a consolidated feature of cosmology Unfortunately, 
the physical mechanism from which this cosmic speed up 
originates is still unclear; the common way to deal with 
this is to assume, in addition to the standard matter 
term, the existence of a further exotic fluid which influ- 
ences the dynamics of the Universe [BHD]. Due to the 
lack of knowledge on the physical nature of this fluid, 
we usually refer to it as dark energy (DE). So far, DE is 
only observationally witnessed, while the micro-physics 
behind it remains totally undisclosed f9lfl3|. One of the 
most dubious properties of DE is that it exhibits a nega- 
tive equation of state (EoS) parameter, counteracting the 
attractive action of gravity [l4[ . The need of a negative 
pressure hints at the non-baryonic nature of DE, since 
no common matter is expected to show such a property. 
Besides that, the total amount of cosmological matter in 
the Universe appears to be dominated by a non-baryonic 
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(cold) dark matter (DM) component, which accounts for 
about 23% of the total energy content of the Universe. 
On the other hand, the common baryonic matter in the 
Universe only accounts for 4% of the whole energy con- 
tent. This shows that the standard visible matter is actu- 
ally not enough to guarantee the stability of structure ob- 
served at different astrophysical and cosmological scales, 
and DM cannot be ignored for the dynamics of the whole 
Universe EGJ]- 

Consequently, our knowledge of the correct cosmologi- 
cal model seems to be lacking of some ingredients. How- 
ever, in order to investigate the effects of DE and DM 
in Einstein's equations, one introduces a common energy 
momentum tensor with a pressureless term, i.e. P m = 0, 
describing the total visible and non-visible matter con- 
tent, and an additional term with a negative pressure 
to represent DE 17J. Together with these assumptions, 



one generally considers a homogeneous and isotropic 
Universe, depicted by the Friedmann-Robertson- Walker 
(FRW) metric, ds 2 = -c 2 dt 2 + a{t) 2 {dr 2 /{\ - kr 2 ) + 
r 2 sin 2 8d(j) 2 + r 2 d0 2 ). In addition, observations of the 
large scale geometry of the Universe suggest a spatially 
quite flat Universe, so hereafter we will assume k = 0. 
To account for the effects of DE, the simplest and most 
tested assumption deals with the introduction of a cos- 
mological constant term A into the Einstein equations. 
According to quantum field theory, the constant is in- 
terpreted as a vacuum energy contribution, and natu- 
rally leads to a negative EoS parameter with a positive 
energy density and negative pressure. The correspond- 
ing model, which is straightforwardly derived by solving 
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the Einstein equations with the cosmological constant, is 
known as ACDM [l8j], a model which by now achieved 
the status of the standard cosmological model. The rea- 
son which induces cosmologists to assume this particular 
model to be the standard one is that it excellently fits 
all observational data with high precision [l9|, [2(| ■ More- 
over, it is only relying on a remarkably small number of 
cosmological parameters, without any ad hoc additional 
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Unexpectedly, observations show that both 
the magnitudes of matter and A are comparable at our 
time. Indicating with Q m and the magnitudes of mat- 
ter and DE respectively, observational bounds show that 
fl\/fl m m 2.7. This feature implies a strange and unex- 
pected coincidence problem -because DE is expected to 
evolve separately from matter, it is quite astonishing to 
imagine that near the present time the two magnitudes 
should be so close to each other. 

On the other hand, another uncomfortable shortcom- 
ing plagues the standard model. The observational limits 
on the magnitude of the cosmological constant disagree 
with the predicted value for about 10 123 orders of magni- 
tude, leading to a serious fine-tuning problem (22|. This 
deeply disturbs the otherwise appealing picture of a cos- 
mological constant and, together with the coincidence 
problem, dramatically afflicts the standard cosmological 
paradigm. Despite its success in explaining the observa- 
tional data, the ACDM model is therefore theoretically 
incomplete 1 , or at least not well understood. Motivated 
by these defects, a mare magnum of different models has 
been proposed during the last decades; as a short sample 
see [25| and references therein. In this work we wonder 
whether the cosmological constant must be considered as 
the real unique explanation of DE, or if there exists a hid- 
den mechanism behind the nature of the cosmic speed up. 
An enticing way to understand if ACDM is the favorite 
candidate for DE is represented by the analyses through 
model independent tests. Such procedures should be able 
to disclose the fundamental nature of DE without postu- 
lating a certain model a priori. In this way, it would be 
possible to analyze the dynamics of the Universe without 
imposing a cosmological constant from the beginning. If 
A really exists, no significant deviations from a constant 
EoS must be found by model-independent tests. As a 
consequence, it is necessary to inquire how much of mod- 
ern cosmology is really independent of the Friedmann 
equations [2a, [28j]. in other words, distinguishing be- 
tween kinematics and dynamics is viewed as a tool to 
discriminate fairly among models, in order to reveal the 
correct cosmological paradigm. 

Surely one of the most powerful model independent 
approaches is represented by cosmography [29]. Cosmog- 
raphy, sometimes also referred to as cosmo-kinetics, was 



1 Moreover, nearly all the extensions of it appear to fail as well. For 
a recent and mentionable alternative, which naturally extends 
ACDM in general relativity, conforming to all the experimental 
bounds, see EE El 113. 



first discussed by Weinberg in [19| and then extended by 
Visser in [26|. The underlying philosophy of cosmogra- 
phy is to involve the cosmological principle only. So, the 
FRW metric is the only ingredient that cosmography uses 
for obtaining bounds on the observable Universe. Cos- 
mography permits us to infer how much DE or alternative 
components are required in regard to satisfy the Einstein 
equations. The idea is to expand some observables such 
as the cosmological distances or the Hubble parameter, 
into power series, and relating cosmological parameters 
directly to these observable quantities. In doing so, it is 
possible to appraise which models behave fairly well and 
which ones should be discarded as a consequence of not 
satisfying the basic demands introduced by cosmography. 
So cosmography strives for the development of a proce- 
dure able to constrain the kinematics of the Universe. 

In this paper, we present an extension of the promis- 
ing approach debated in Visser et al. [29j ; we devote 
our efforts both to constraining ACDM and investigating 
whether it is the only possibility to explain the cosmolog- 
ical acceleration, or whether there are prominent alter- 
natives [30j . In particular, we adopt the idea of cosmog- 
raphy developed in [26l [29l [3ll - 434| and we improve it, by 
assuming an extended class of fitting quantities, assem- 
bled by a number of different cosmological distances. For 
theoretical reasons, which we will discuss in the next sec- 
tions, we introduce new parameterizations of the redshift 
variable besides z, in order to improve the fitting proce- 
dure. These parameterizations are designed to reduce 
the problems associated to the experimental analysis at 
redshift z > 1. Thence, we make use of the most recent 
data of the Union 2 supernovae la (SNela), of the Hubble 
Space Telescope (HST) measurements of the Hubble fac- 
tor, and of the H(z) compilations |3]|, through a Markov 
Chain Monte Carlo (MCMC) method, by modifying the 
publicly available code CosmoMC [36|. Afterwards, we 
also include a parametrization of the cosmological dis- 
tances in terms of the EoS of the Universe as a whole and 
of their pressure derivatives. This allows us to directly 
fit the EoS of the Universe without having to undergo 
disadvantageous error propagation. In doing this, we as- 
sume the validity of the cosmological principle, and of 
General Relativity, since for this analysis it is necessary 
to invoke the Friedmann equations. This gives us certain 
constraints on the EoS and on the pressure derivatives in 
the framework of General Relativity. 

The paper is organized as follows. In Sec. II we dis- 
cuss the physics behind the concept of cosmography and 
its implications for modern cosmology in more detail; we 
introduce the new cosmographic coefficient m, and dis- 
cuss how to build up a viable alternative parametrization 
to the redshift z. We then propose three new parame- 
terizations and we study their properties, in view of the 
fitting procedure. In Sec. Ill we apply these recipes and 
present the results of the cosmographic fits by a numeri- 
cal MCMC analysis. Section IV deals with the concept of 
the EoS; in particular we relate the EoS and the deriva- 
tives of pressure to the luminosity distance, in order to 



3 



carry out a direct fit to cosmological data. We use these 
results to derive constraints on ACDM and the most 
generic DE model, characterized by the function G(z), 
reducing to G(z = 0) = 1 — O m at z = 0. In Sec. V fi- 
nally we draw conclusions and give an outlook to further 
paths of investigations. 



II. THE ROLE OF COSMOGRAPHY 

In this section we focus on the role of cosmography in 
modern cosmology. Its aim is the study of the kinematic 
quantities, characterizing the cosmological scenario. For 
this reason, cosmography is also called cosmo -kinetics, 
or kinematics of the Universe. We therefore limit our- 
selves to the smallest number of assumptions possible. 
First, we presume the validity of the cosmological prin- 
ciple. Second, we suppose that the EoS of the Universe 
is determined by a non-specihcd number of different cos- 
mological fluids 2 . We assume that the total pressure of 
these fluids, namely P, can be written as P = ^2i;P%, 
and its total EoS parameter uj = J2i Pi/ J2i Pi- Here 
the index i runs over all the involved cosmological fluids. 
Following these assumptions, the par adig m of cosmogra- 
phy was first developed by Weinberg Qjl , who proposed 
to expand the scale factor in terms of a Taylor series 
around the present time to 3 . Following these recipes, 
it is natural to expect that many other physical quanti- 
ties of interest, apart from a(t), can be expanded as well. 
The power series coefficients in the expansion of the scale 
factor are known in the literature as cosmographic series 
(CS), when evaluated at our time fo; these quantities are 
related to the scale factor derivatives. 

A feasible consequence of this prescription is that cos- 
mography does not depend on the choice of a cosmolog- 
ical model. As a matter of fact, almost all cosmological 
tests assume a priori that the model under considera- 
tion is statistically favored; however, this creates a de- 
generacy among models and often it remains difficult to 
understand which model is really favored. Cosmography 
is, among various cosmological tests, one of the ways to 
alleviate that degeneracy. However, although cosmogra- 
phy is reviewed as a model independent procedure, a few 
words should be spent regarding the role of the spatial 
curvature, k. In particular, modern cosmological data are 
not enough at present to fix stringent convergence limits 
on the CS and k. It is possible to show that the term 
proportional to the (present) "variation of acceleration" , 



These fluids include matter, radiation, curvature, dark energy 
and so forth. 

Instead of the scale factor, also the Hubble parameter or the 
luminosity distance could be expanded. At late times it is al- 
lowed to neglect radiation within the energy momentum tensor 
in Einstein's equations. In addition, as will be clarified later, 
we superimpose a spatially flat geometry, in accordance with the 
most recent observations. 



i.e. jo, cannot be measured alone. Defining the curva- 
ture density as fi^ = Q — 1, where Qq represents the 
total density of the Universe, then one measures jo + Qo 
13. Motivated by WMAP 7 results we propose 

here to restrict the analysis to the spatially flat case, in 
which k = 0. This naturally overcomes the dependence 
of jo on r2o, letting cosmography be independent of any 
particular cosmological framework. 

Now we have all the ingredients to expand the scale 
factor into a series, yielding 
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where we truncated the series at the sixth order in 
At = t — to. Here, we assume that t — to > 0. More- 
over, the constant ao is the scale factor evaluated today. 
Without loss of generality, it is licit to identify hereafter 
ao = 1- Equation (JXJ) can be recast as 
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with the definition of the cosmographic coefficients as 
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Having an expansion for a(t) is equivalent to having an 
expansion of the redshift z, in terms of H At (see [38|). 

As previously stressed, Eqs. (O, if evaluated at our 
time, are referred to as the CS. The subscript "0" in Eq. 
@ indicates that the coefficients are evaluated at t = to- 
In particular, each term has its own specific physical in- 
terpretation. For example, q, the so-called acceleration 
parameter, specifies whether the Universe is accelerating 
or decelerating, depending on the sign. An accelerat- 
ing Universe leads to —1 < go < 0. On the contrary, a 
positive jo implies that g changes sign as the Universe 
expands, and so forth for all the rest of the parameters. 
We usually attribute the names of jerk and snap to j and 
s respectively; so far, no universal name is associated to 
I. Here, we additionally introduce to as a further higher 
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order term. It is useful to combine the CS among them- 
selves and express them in terms of each other, yielding 
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3 [21 + 5 (24 9 + 18g 2 + 2q 3 - 2s - qs + 8)] . 

Here, the dots and the numbers in brackets indicate 
the derivatives with respect to the cosmic time. By con- 
verting the derivatives in Eq. Q from time to redshift, 
and inverting the relations, we can obtain the Hubble pa- 
rameter as an expansion in terms of redshift, H(z) (the 
results can be found in Appendix[B]). Analogously, we can 
also expand other observable physical quantities in order 
to fit the cosmological data with the obtained functions. 

It would be interesting to expand commonly used no- 
tions of cosmological distances to the same order of the 
Taylor expansion as the scale factor before. To this end, 
let us now introduce several examples of distances be- 
tween two objects in cosmology, following the prescrip- 
tions given in [31j. Here, we state the luminosity distance 
cIl and other four alternative distances, namely the pho- 
ton flux distance dp, the photon count distance dp, the 
deceleration distance d,Q and the angular diameter dis- 
tance <1a- These distances arc defined as 



d L 
dp- 
dp 

d Q 
dA 



a r (l + z) = r • — — 
a(t) 



d L 



1 



ro 



(1 + 2) 1/2 " y/a{t) 



d L 



(5) 



(1 + z)3/ 2 

d L 



(1 + z) 2 



= r ■ a(t) 



The last four notions of distances are less commonly 
used in literature. Besides the luminosity distance dL, 
which gives the ratio of the apparent and the absolute 
luminosity of an astrophysical object, we consider the 
photon flux distance dp, which is not calculated from 
the energy flux in the detector, but from the photon flux, 
which is experimentally easier to measure. The photon 
count distance dp is based on the total number of pho- 
tons arriving at the detector as opposed to the photon 



rate. The so-called deceleration distance dQ has been 
introduced in [26j without having an immediate physi- 
cal meaning, but in return a very simple and practical 
dependence on the deceleration parameter qq. Finally, 
the angular diameter distance dA was defined in [l^] as 
the ratio of the physical size of the object at the time 
of light emission and its angular diameter observed to- 
day. To completely determine the distance expansions, 
we still need to calculate ro. It is defined as the distance 
r a photon travels from a light source at r — ro to our 
position at r = 0. It is defined as 4 



ro 



dt' 

W) 



(6) 



We can calculate this quantity by inserting the power 
series expansion for the inverse of the scale factor and 
integrating each term in the sum separately. Finally the 
results are used to complete the calculations of the cos- 
mological distances in terms of the redshift z. We report 
in Appendix IA II the expansions of all the distances in 
Eqs. ((SJ) in terms of z. These results can be compared 
with those of [29 ], in which the authors truncated the 
series at a lower order. In particular, they claimed the 
need of using all the distances for a cosmographic test. In 
principle, this may be true, because all the various cos- 
mological distances rely on the fundamental assumption 
that the total number of photons is conserved on cosmic 
scales. Hence, there is no reason to discard one distance 
for another one, since all of them fulfill this condition. 
Unfortunately, there exists a duality problem plaguing 
such distances [HI . This problem is so far an open ques- 
tion of observational cosmology |4(J. On the other hand, 
it has been suggested that the luminosity distance d^ is 
fairly well adapted to the cosmological data used in com- 
bined tests with supernovae la and HS T l4l| - |43j . Even 
though this topic is still object of debate |44{. unlike Cat- 
toen and Visser [2i| we limit our attention to dL only. We 
motivate this choice with the above considerations on the 
good adaptation of dj, to the data, and with its general 
use in literature 
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There are two main problems arising in the context of 
cosmography. In principle, the Taylor series is expected 
to diverge at z > 1. This is a consequence of the fact 
that we are expanding around z ~ and so when z > 1, 
we get problems with convergence. Moreover, the finite 
truncations we made represent only an approximation of 
the exact function, giving therefore possibly misleading 
results. Thus, while the second problem can be alleviated 
by expanding to higher orders in adding more coefficients, 
this measure can introduce divergences into the analysis. 
These issues are intimately connected to the problem of 



4 Here we have omitted a factor of c in the numerator; for now and 
for the rest of the theoretical calculations in this paper, we will 
assume c = 1. 
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systematic errors. In fact, if errors are large enough, it is 
possible that bad convergence may afflict the numerical 
results. 

We improve the accuracy of our work by using the 
Union 2 compilation, which reduces the problem of sys- 
tematics, easing the second problem. On the other hand, 
in order to overcome the first issue, different parameteri- 
zations of the fitting functions can be taken into account. 
The idea is to carry out the expansion with a different 
variable which is used ad interim and is constructed to 
be limited in a more stringent interval. While z G [0, oo], 
a new variable should for example be restricted to the 
interval [0,1]. 

A fairly well-known possibility is represented by the 
variable 



2/i = 



1 



(7) 



frequently used in literature [2 a . 1271 ] . The limits in the 
past Universe, i.e. z G [0, oo], read y\ G [0,1], while in 
the future, i.e. z G [—1,0]: y\ G [— oo,0]. Immediately 
we notice that y\ can be expanded as z before as y\ = 
yi(HoAt); then, it is feasible to invert it, having H^At 
in terms of y\. Then, we can express the distances as 
functions of 2/1 (for the results, see Appendix I A 2|) . 

Furthermore, we propose a way to construct other vi- 
able parameterizations of the redshift variable. To this 
end, we introduce below three new propositions, namely 
2/2,2/3 and y 4 , as 



j/2 = arctan 
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whose limits are, for z G [0, oo] : y% G [0, f],j/3 G 
[0,0], 2/4 G [0,f] and z G [-1,0] : j/2 G [f ,0],t/ 3 6 
[— i,0],2/4 G [— f ,0] and in which we used the definition 
of the scale factor, i.e. a = (1 + z) -1 . 

We adopted the arctan in the parameterizations of 2/2,4 
because it behaves smoothly and it is suited to give well- 
defined limits at z — » 00. On the contrary, 2/3 is a poly- 
nomial in z; so apparently, we would not expect it to 
lead to significantly different fitting behavior, but just to 
represent an alternative worth investigating. 
Equations (JSJ can be expanded into Taylor series for 
z < 1, and then be inverted for 2/2,3.4 in order to give 
an expression for i/oA/;(2/2.3,4), up to sixth order. Then 
the distances as functions of 2/1,2.3,4 can be calculated as 
well (results see Appendices IA 3tiA 5|) . 

By definition, all these parameterizations are built up 
to avoid divergences at z > 1. Thus, one can wonder 
whether all of them turn out to be equally suitable for 
constraining the CS. The answer can be partly predicted 
by comparing the supernova data of the luminosity dis- 
tance, as in Fig.[TJ for z and 2/1,2,3,4- The worst example 
is clearly the redshift 2/3- Its definition suggests that it 



scales down more quickly, compared to the redshift z, 
as can also be seen in Fig. [TJ This means that a re- 
gion of z G [0, 1.5] is reduced to a much smaller interval 
2/3 G [0,0.5]. Thereafter, we expect that, when the curve 
bends too quickly, the fits become more difficult. It fol- 
lows that, as the curve trends become more extreme, a 
suppression of lower redshifts to the advantage of higher 
ones can occur. In other words, z G [0,0.5] weighs less 
than z > 0.5; therefore, we guess that 2/3 would work 
better if all the cosmological data were for z> 1. 

As it can be seen in Fig. [1] the luminosity distance 
curves of data points over redshift are slightly flexed, be- 
coming steeper towards higher redshifts. Also the red- 
shifts 2/1, 2/2 and 2/4 lead to steeper curves than z, how- 
ever, redshift 2/3 behaves the most extreme. According 
to the cited criteria, 2/3 is the least suitable of redshift 
notions. This conjecture is also backed by the fitting re- 
sults, which confirm that 2/3 does not work well in the 
application to SNela data. Another disadvantage of 2/3 
is that it docs not have a uniquely defined inverse. For 
these reasons, we decided to take it out of the analysis. 
Through similar arguments, we remove the second- worst 
redshift, 2/2, as well. Summing up, in order to contrive a 
viable redshift parametrization, the following conditions 
must be satisfied: 

1. The luminosity distance curve should not behave 
too steeply in the interval z < 1. 

2. The luminosity distance curve should not exhibit 
sudden flexes. 

3. The curve should be one-to-one invertiblc. 

From Fig. [T] we notice that the last introduced redshift 
2/4, although still producing a steeper curve than z, is ex- 
pected to work better than 2/1 . Thus, the rest of the anal- 
ysis, including further calculations and fittings, is carried 
out for the redshifts z, 2/1 and 2/4 ■ 



III. THE FITTING PROCEDURE AND THE 
COSMOGRAPHIC RESULTS 

In Sec. II, we explained how to construct viable cosmo- 
logical parameterizations to alleviate the problems asso- 
ciated to cosmography. To this end, we introduced 2/4 and 
we investigated its theoretical viability to fit the cosmo- 
logical data. We now have all the ingredients to develop 
a MCMC procedure to find cosmographic constraints for 
the CS, using the three redshifts z, 2/1 and 2/4 for the fit- 
ting analysis. We will include the sixth order of the CS, 
i.e. mo, and particularly focus on the following aspects: 

• We investigate whether redshift 2/4 is actually suit- 
able to obtain accurate values of the CS, as the- 
oretically predicted, and we explore the ranges of 
low and higher redshifts; 
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redshift 



FIG. 1: (color online) Luminosity distance (in units of 
10 26 m), over different redshifts z (red), y\ (green), yi (or- 
ange), J/3 (blue) and 3/4 (black). 



• We analyze how well the introduced CS parameter 
?no can be constrained, in particular, we examine 
whether its introduction significantly enlarges the 
dispersion of the estimation of the other parame- 
ters; 

• We find out if the concordance model, i.e. ACDM, 
is in agreement with the cosmographically found 
numerical results. 

For our purposes, we use the data of the SNela Union 
2 compilation by the supernovae cosmology project [42j]. 
We also adopt the HST measurements on 600 Cepheides, 
which impose a Gaussian prior on the Hubble parameter 
today of H = 74.0 ± 3.6km/s/Mpc 0, and the mea- 
surements of the Hubble parameter H(z) at twelve dif- 
ferent redshifts ranging from z = 0.1 to z = 1.75 [47} . We 
divide our analysis into two sets of observations, namely 
set 1, which comprises Union 2 together with HST, and 
set 2, being the Union 2 dataset with both HST and 
H(z) measurements. For the sake of completeness, it 
is in order to cite recent works using gamma ray bursts 
(GRBs) as possible distance indicators [HI, |33|, [48| . How- 
ever, considering GRBs in such analyses is mere specu- 
lation, since GRBs are not standard candles j4§|. This 
approach seems to result in wrong estimations, or at least 
inadequate results. It will be shown that our results dif- 
fer from those obtained by using GRBs; we will show a 
set of results in better agreement with ACDM than those 
obtained by using GRBs, which points to an inadequacy 
of the use of GRBs in cosmography. We also exclude ob- 
servation data from baryonic acoustic oscillations from 
our analysis; we deem that introducing baryon acoustic 
oscillations means reducing the model-independence of 
the whole analysis [50j | . 



In the following, we will use the CS combined together 
in three sets with different maximum order of parameters: 



A = {#0,90, jo, So}, 

B = {H ,q ,jo,s ,lo} , (9) 

C = {H Q ,q ,j ,s J ,m }. 

We expect a slower convergence of the last dataset, 
since the introduction of mo can decrease the accuracy 
of convergence. To constrain the parameters, we use a 
Bayesian technique in which the best fits are those max- 
imizing the likelihood function C oc cxp(— x 2 /2). Since 
the different observations are not correlated, the function 
X 2 is simply given by the sum \ 2 = Xunkm 2 + Xhst + 
Xh( z ) ■ We explore the space of parameters with a MCMC 
approach, modifying the publicly available code Cos- 
moMC !36]. We do the analysis for the three sets' param- 
eter space, for the two sets of observations, and for each 
of the three considered redshifts. Accordingly, we per- 
form 18 different constraint parameter analyses. To ob- 
tain the posterior samples we assume flat priors over the 
intervals -6 < q < 6, -20 < j < 20, -200 < s < 200, 
-500 <l < 500, and -3000 < m < 3000. 

In Tables HI HU and IIIII we show the best fits and their 
ler-likclihoods for the redshifts z, y\ and 7/4 respectively 

In Fig. [5] we compare the 1-dimensional marginalized 
posterior distributions for each parameter and each red- 
shift for set 2 of observations. We note that the param- 
eters Iq and mo are not well constrained when using the 
redshift y±. By introducing 7/4 it is possible to overcome 
this issue, obtaining good results on Iq and mo as well. 
As expected, the redshift z appears to be statistically 
more favored than j/4. We confirm what we conjectured 
in Sec. II: the sixth order of the Hubble expansion in z 
works better than any other paramctrization, if z < 1. 
Nonetheless, 7/4 should be taken seriously as a possible 
alternative to y±, when z > 1. Similar conclusions have 
been drawn for the set 1 ; see Tables HI UH and IIIII The 
only caveat is that for set 1 the results appear to be less 
accurate as in the previous case. 

We proceed to determine if the introduction of the pa- 
rameters Iq and mo is convenient. To this end, in Fig. [3] 
we plot the first four CS parameters' posterior distribu- 
tions for the three parameter sets, comparing their statis- 
tical widths. We note that the dispersions are enlarged 
considerably when we add the ^o parameter; however, 
introducing mo does not substantially broaden the pos- 
terior distributions; besides, the standard deviations of 
the posteriors are in a proportion 1 : 2.28 : 1.92 for jo 
and 1 : 5.66 : 8.38 for sq. 

In Fig. 2] we present the summary of the results for 
the redshift 7/4 and parameter set C by plotting the 2- 
dimcnsional contours and the likelihood samples. 
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TABLE I: Table of best fits and their likelihoods (la) for redshift z, for the three sets of parameters A, B and C. Set 1 of 
observations is Union 2 + HST. Set 2 of observations is Union 2 + HST + H(z). 



Parameter 


A, Set 1 


A, Set 2 


B, Set 1 


B, Set 2 


C, Set 1 


C, Set 2 




Xmin = 530.1 


545.6 


530.1 


544.5 


530.0 


544.3 


Ho 

qo 

jo 
s 
lo 

III,, 


74 35+ 7 - 30 
-0 70 85+ - 6074 

U. 1 WJ_ .50.52 

1.605+2 

— 4.481 
9 5^ + 60. 61 


74.22+i;2i 
-0.6149+° St 
1.030+°- 7 ^ 

— 1 .001 

16+ 1 - 45 

u ' iu -1.03 


73.77±|:1| 
-0.6250+° :f 9 l° 
0.392+*^5 

— 4.511 
— 5 5Q+33.74 

° ,dU -89.19 


74.20+.!:°! 
-0.6361tgl-g 
0.994+1 ™* 

— 2.665 

— 1 47+4.20 
4 47+41.53 


73.72+ 7 - 47 
-0.6208+°: 4 7 « 
-1.083+1^9 

— 2.218 

-25.52+«;Sg 
N.C. 
N.C. 


73.65±!:f 6 

-0.5856t8; 3 445 
-0.117+3-^1 

— 1.257 

— 7 71 +14-77 
' ' ' -7.83 

C CI^ + 23.39 
°- d0 -27.86 

71 qS + 382.17 




b. N.C. 


means the results are not cor 


Notes. 
ff is given in Km/s 
lclusive. The data do 


/Mpc. 

not constrain the parameters 


sufficiently. 




TABLE II: Table of best fits and their likelihoods (1 
observations is Union 2 + HST. Set 2 of observations 


j) for redshift j/i, for the three sets of parameters A, B and C. Set 1 of 
is Union 2 + HST + H(z). 


Parameter 


A, Set 1 


A, Set 2 


B, Set 1 


B, Set 2 


C, Set 1 


C, Set 2 




Xmin ~ 530.1 


550.1 


529.9 


544.5 


530.0 


545.1 


Ho 

qo 

Jo 
so 
lo 


74 05+ 7 - 90 
-0.6633i°-= 7 f o 
1.268±|if| 
1.21+|^ 


75.25±1;^ 7 
-0.4106+8;"" 
-7.746+J"* 
-88.91±K:g 


73.68+. 7 - 77 4 

-o.ooo4+5:^i3 

-13.695+f^ 1 
-180.95+^ 
N.C. 


73.30+i-g 
-0.2652+° 07 ° 77 
-7.959+£tf 
-112.63+^ 6 f 
jV.C. 


73.91± 7 :™ 
-0.536018;I 4 9 6 6 8 5 
-1.646+^ 3 f/ 
-30.97±1°: 47 
N.C. 


74 4Q+5-07 
-0.4624+° o f« 
-1.862+ii 3 ° 97 i 

— 16 95+ 73 68 

±u.»j_ 3g 7g 

N.C. 


mo 










N.C. 


N.C. 



Notes. 

a. /f n is given in Km/s/Mpc. 
b. N.C. means the results arc not conclusive. The data do not constrain the parameters sufficiently. 



IV. THE CONNECTION BETWEEN THE CS 
AND THE EOS OF THE UNIVERSE 

In Sec. I, we outlined that an expression for the EoS 
is naturally associated to each cosmological fluid in a 
given cosmological model; in particular, in thermody- 
namics the EoS characterizes the properties of such a 
fluid. Under the assumption of a specific gravitational 
model, the quest of understanding the expansion history 
of the Universe is equivalent to reproducing the correct 
EoS during different phases. Many physical mechanisms 
are hidden in the EoS parameter lo. Finding the correct 
EoS has thus high importance in cosmology, because it of- 
fers a key to understanding the micro-physics associated 
to DE and/or DM. In this section we will assume that 
General Relativity gives a correct description of gravity 
at the scales under consideration. 

Under the hypotheses of cosmography, we cannot a 
priori assume an EoS of the Universe, because we were 
not specifying any particular cosmological model at the 
beginning. We recall that the EoS of the Universe is 
given by ui = Pi/ pi, where the subindcx i refers to 
the different fluids that the Universe comprises. Hence, 



to evaluate lj, one needs to know the total pressure, 
P = ^ i Pi 1 and the total density, p = X^/ 9 *- Even 
though we do not assume any EoS explicitly, it is possi- 
ble to expand the pressure in terms of the cosmic time or 
redshift variables, i.e. z,y\,y&. By expanding the pres- 
sure into a series, it is possible to predict the values of 
its derivatives with respect to the cosmic time or the red- 
shift variables. In fact, one can relate the derivatives of 
P to the CS; therefore, by substituting the values of CS 
in terms of these derivatives into the luminosity distances 
we are able to directly fit the parameters of the EoS of 
the Universe from luminosity distance data. 

The reason of constraining the pressure derivatives lies 
in the possibility of discriminating among models; in 
principle, a model which does not satisfy such bounds 
can be easily discarded. The expansion of P in terms of 
the cosmic time is formally given by 

where yi — z, y±, j/4. By truncating the series at the 
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TABLE III: Table of best fits and their likelihoods (la) for redshift j/4, for the three sets of parameters A, B and C. Set 1 of 
observations is Union 2 + HST. Set 2 of observations is Union 2 + HST + H(z). 



Parameter 


A, Set 1 


A, Set 2 


B, Set 1 


B, Set 2 


C, Set 1 


C, Set 2 




Xmin = 530.3 


544.8 


529.7 


544.6 


529.9 


544.5 


H 


74.55+?;K 


73.71+11? 


73.95+™ 


73.43+1;?! 


74.12+.?; 2 ? 


73.27+S;«" 


'in 


-0.7492^0.5899 


-0.6504+2-4275 


-0.461 1+0;«» 


-0.7230+8;»gi 


-0 4842+ 2 - 7126 


-0.7284±g : S°|2 


jo 


2.558+^ 




-3.381+J»ig» 


2-017±l;i« 


-1.940+1;?« 


2 148+ 3 - 414 

"■ A ™— 4.038 


s 


O atC+74.69 
a - aa -26.69 


3.151+?;™ 


—37 67+ 89 - 51 

J ' ,u ' -60.10 


5.278+ii : og 


-13.48±g;SS 


9 1 7Q+42.126 
^•- L 1 ^-35.919 


lo 






AT.C. 


—0 1 S+ 96 - 75 


N.C. 


— 11 60+ 193 SS 


m 










N.C. 


70 + 2497.8 
' U ' -2254.5 




b. N.C. n 


leans the results are not co 


Notes. 

K is given in Km/s/Mpc 
lclusivc. The data do not c 


jnstrain the paran 


otcrs sufficiently. 





fourth order, and the continuity equation 
dp 



dt 



3H(P + p) = 0, 



(11) 



and the Friedmann equation H 2 = ip, we can write down 
an explicit dependence of the pressure and the coefficients 
on the CS as follows 5 



In addition, by using Eq. ([TTj) . Eqs. (fl"2")) and the Fried- 
mann equation H 2 = ^p, we find the expression for the 
EoS parameter of the Universe as 



1 



(14) 



A. Fitting the EoS 



P = -H 2 (2q - 1) , 



dP 

H 
d 2 p 

~dt 2 
d 3 P 
~dt? 

d 4 P 



\h 4 (j -3q-s-3), 



(2s +j-l + q(21- j) + 6q 2 + 12 
j 2 +31 -m- U4q - 81q 2 



3 

6q 3 - 12j (2 + q)- 3s - 3qs - 60 



(12a) 
(12b) 
(12c) 

(12d) 
(12c) 
(12f) 



where we evaluated the derivatives up to the order of ttiq. 
We list the coefficients 4-f, 4~f and 4~f in Appendix ICl 

az dy 1 ay 4 

For completeness, we write down the transformation laws 
between z, 2/1,1/4 and the cosmic time t, i.e. 



d_ 

at 



[l + z)H- 



d_ 

dz 
d 

dyi 



(13) 



cos 2/4 (cos j/4 + sin 2/4 ) H 



5 Here in Eqs. (1 1 L'[) we are considering = c = 1 f or brevity. 

These factors are considered again in the numerical simulations. 



In this subsection, our goal is to obtain constraints on 
the EoS and the pressure derivatives by inverting Eqs. 
(|12b[) - (|12fj) and by rewriting the luminosity distance as 



a function of ojq, 3^7 

In other words, we use Eqs. 
cated henceforth by 



dy'f 

pi) 



d 4 p 

dy 4 

flU 



and Ull), indi- 



V = 



ay; 



2^3 



3 ' P 4 



d 4 P 



to express the vector CS = {qo, jo, so, !o,mo} as a func- 
tion of the EoS parameter and the pressure derivatives, 
CS = CS(ujq, Pi, P2, P3, Pa). The purpose is to plug 
those results into the expressions for the luminosity dis- 
tance in order to find numerical best fit values of the 
new set of parameters, using d^ = di(wo, Pi, Pi, P3, P4) 
in the numerical analysis. 

In principle, there exists also an alternative procedure, 
which consists in taking the results already obtained for 
the CS and to propagate the errors through Eqs. (|12b[) - 
(|12fp , without performing another fitting procedure. But 
in choosing this way, we would face an unacceptable 
increase in the errors, as opposed to the direct fit of 
{too, Pi, Pi, P3, Pi}. Thus, in order to reduce the error 
propagation, the simplest and most straightforward way 
is to evaluate the coefficients by a direct fit of the lumi- 
nosity distance. The search for the best-fit values for the 
new set of parameters is performed by using the proce- 
dure of MCMC simulations developed in Sec. III. 

For statistical reasons, we choose set 2 of observations, 
which is more complete and suitable for this kind of fit. 



64 66 68 70 72 74 76 78 80 82 
Ho 









i Model A 

Model li ■ 
V Model C 


.-.-Li- i 





-60 -40 -20 20 40 60 




FIG. 3: (color online) 1-dimensional marginalized posteriors 
for Ho, go, Jo and so, using set 2 of observations (Union 2 
+ HST + H(z)). Solid (red) line is parameter set A, dotted 
(green) line is parameter set B and dashed (black) line is 
parameter set C. 



The explicit expressions for the luminosity distance in 
terms of the different rcdshift parameters are reported 
for completeness in Appendix [Dj As in Sec. Ill we limit 
our analysis to d^- 



FIG. 2: (color online) 1-dimensional marginalized posteriors 
for the complete CS (parameter set C), using set 2 of observa- 
tions (Union 2 + HST + H(z)). Dotted (green) line is redshift 
z, dashed (black) line is y\ and solid (red) line is 3/4. 



Figure [5] shows the obtained marginalized posteriors 
and in Tab. IIVI we present the summary of the results. 
We observe the same hierarchy of rcdshifts as in analysis 
of Sec. Ill; meaning that our new "redshift" introduction, 
?/4, appears to be statistically favored with respect to y\. 
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TABLE IV: Table of mean values of the posteriors and their likelihoods (la) for the three redshifts, using set 2 of observations 
(Union 2 + HST + H(z)). Pi(z), Pi(yi), Pi(y 4 ), fbO), Pi(.Vi) and P 3 (z) are in units of 10 4 c 2 /k, P 2 (yi) and P 4 (2/ 4 ) in units 
of 10 5 c 2 /k, P3(yt) and Pa(z) in units of 10 6 c 2 /k, Pj(yi) in units of 10 7 and Pi(yi) in units of 10 8 c 2 /k. 



Parameter 


Redshift z 


Redshift yi 


Redshift y 4 


H 

Pi 
Pi 

Pz 

Pi 


74 23+ 2 31 

1 -2.36 

— 71 74+0.0922 
u ' 1 1 -0.0964 

—0 20Q+o- 347 
0.630+t°« 

-0. 10710.090 


74.20±m 
-0.7439+°;^ 
-0.991±3;lfl 
-0.134+i;«g 
0.205±° o: ^ 
—0 150+ - 209 


75.70± 2 2 ;f 6 
-0.7315±g;li~ 
-0.228+° ; f 2 i 
-0 246+ 4133 

u "" u -3.927 

9 1 7+0.625 
u -^r ' -0.400 

-0.289+S : ;g 



Notes. 

a. H Q is given in Km/s/Mpc. 



-0.4 
-0.6 
-0.8 



1 

70 75 81 

# 




-20 2040 -200 200 
2000 { 



-2000 1 | -2000 

70 75 80 -1 .20.804 



-20 2040 -200 200 -2000 2000 



FIG. 4: Marginalized posterior constraints for redshift y 4 and 
parameter set C, using set 2 of observations (Union 2 + HST 
+ H(z)). The shaded region and the dotted lines show the 
likelihoods of the samples. 



B. Comparison with models 

The results found by fitting the EoS show no conclu- 
sive evidence for a pressure varying in time, as shown by 
the la confidence levels in Tab. IIVI This means that a 
negative constant pressure model is favored for depicting 
the cosmic speed up. Thus, the only two models account- 
ing our results appear to be the concordance model and 
the vanishing speed of sound model (VSSM), proposed 
in [24|, [3(J . Their EoS parameters (neglecting radiation 



components at late times) are w = —1/(1 + ft m /^Aa 3 ) 
and oj = -1/(1 - (f - fl m /n x )a- 3 ), for ACDM and 
VSSM respectively. They imply wq ~ —0.73 nowadays, 
by using the values of Vt m ~ 0.27 and fix — 0.78, 
£ ~ —0.025. These results have been confirmed by the 
present analysis. 

Using set 2 of observations we worked out the ACDM 
model, which is for our purposes and the redshifts 
involved sufficiently described by the two parameters 
{H, m h 2 ,8}. The philosophy is to estimate Vl m h 2 and 
9 by the Monte Carlo simulation and then substi- 
tute it into the CS in the ACDM model. To evalu- 
ate the CS for ACDM we use Eqs. (0} with H = 
H yJVL m (l + z) 3 + 1 - n m , yielding 



jo 
so 

k 

mo 



-1 
1. 

1 - 

1 
1 



n, 



2 
m 

27 



(15) 



27 



810 



m * 

2 
m 



81 



ft? 



Any significative tension between the CS values ob- 
tained in this way and the values derived before using 
cosmography would be an indication of the validity of 
a different theory, other than the concordance model. 
Table [V] gives the summary of the likelihoods of the 
estimated and derived parameters. Comparing these 
values with the ones obtained for models A, B and C we 
note that all our results are compatible with the ACDM 
model within the limits of error. 

Now, we want to factorize the effects of DE by assum- 
ing a cold dark matter model, for which the Friedmann 
equation is given by 

H = H ^/Q m (l + z) 3 + G(z), 
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-300000 
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FIG. 5: (color online) 1-dimensional marginalized posteriors for the complete set of parameters of the EoS analysis, using set 2 
of observations (Union 2 + HST + H(z)). Dotted (blue) lines are used for z, dashed (black) lines for yi, and solid (red) lines 
are for redshift j/4. 



where G(z) models the corresponding DE term. We want 
to obtain constraints on G(z) and its derivatives with re- 
spect to rcdshifts for the present time z = 0. To this 
purpose we calculate the derivatives of the Hubble pa- 
rameter as stated above with respect to the redshift and 
equate the results to the derivatives of the Hubble param- 
eter in terms of the CS. Thus we obtained the derivatives 
of G(z) with respect to z (results see Appendix [E| ■ The 
same has been done for the redshifts y\ and 2/4. 

We can now use the numerically obtained values for 
the CS and the value Cl m = 0.274^;^ obtained by the 
WMAP7 collaboration [2(3] to calculate the G (k) (y l ) and, 
by error propagation, their la error bars. The results can 
be found in Table IVII With these values it is possible to 
investigate the compatibility of a model of dark energy 
characterized by the function G(z) with the observational 
data at present time. The value of G{z) evaluated today 
can be estimated by the flat space condition, which im- 
plies Go = G(z = 0) = 1 — TO or, assuming Gaussian 
distributions, G = 0.726±g;gi|. 

From Table fVTl we conclude that our results are con- 
sistent with a constant function G(z) = which is the 
case of the ACDM model, or the case of an emergent 
constant as shown in 0, [30[ . 

V. CONCLUSIONS 

In this work we emphasized the importance of con- 
straining the Universe dynamics through the use of a 
model-independent procedure, which does not a priori 
assume the validity of a particular cosmological model. 
In other words, we used the so called cosmography, some- 
times also referred to as cosmokinctics, to investigate the 



TABLE V: Table of best fits and their likelihoods (la) for the 
estimated (top panel) and derived (lower panel) parameters 
for the ACDM model, using set 2 of observations (union 2 + 
HST + H(z)). 



Parameter 


Best Fit (la) 


^ m. h 

e 


1447+0.0181 

-0.0174 

1 060 +OO2 ° 


H 


74.05+™° 


qo 




jo 


1 


so 
lo 


-0.2061+° 1772 

-"--0.2015 
9 7 74+0.485 


mo 


-8.827±l : |Sf 



Notes. 

Hq is given in Km/s/Mpc. Here h is defined through the relation 
Hq = 100 hkm/s/Mpc, and is the ratio of the sound horizon to the angular 



kinematics of the Universe. By following [29(, we per- 
formed an analysis combining theoretical derivations of 
cosmological distances and numerical data fitting, using 
the Union 2 compilation, together with the HST and the 
H(z) samples. For the fitting, we introduced new pa- 
rameterizations in addition to the conventional redshift 
z to express the cosmological distances, underlining their 
importance to avoid divergences at high redshifts and to 
increase the accuracy of the analysis. Moreover, we pro- 
posed prescriptions to build up new viable paramctriza- 
tions, able to overcome these issues. 

We considered three further parameterizations, pre- 
dicting that only one is really a viable option for fitting. 



12 



TABLE VI: Table of derived values and their likelihoods (la) for the derivatives of G(yi) for the three redshifts 2,2/1,2/4, 
evaluated at t = to, using set C of parameters and set 2 of observations (union 2 + HST + H(z)). 



Parameter 



Best Fit for z (la) 



Best Fit for 2/1 (Iff) 



Best Fit for 2/4 (Iff) 



(i) 



(2) 



Go 
G\ 



(3) 



G, 



(4) 


(3) 



0.0068+^f 
-2.22+JII 
13.65±|° 4 

— 17 24+ 35 - 27 

' '"-31. 89 
— 129 74+174.29 



o.25±i;»; 

_4 71+9.13 



7 4 76+ 55 73 



iV.C*. 
iV.C. 



Notes. 
The data 



-0.2788±™ 

1 7^S4+4.SSS27 
J-. 1 t>UM_ 4 g4215 

-3.43039+&JSS 
-16.2906t|™ 
-8.74593±|l|;S?i 



neters sufficiently. 



By using the luminosity distance, i.e. c?l, in terms of the 
standard redshift z, of the alternative parametrization 
yi = and of the newly introduced j/4 = arctan(z), we 
obtained bounds on the cosmographic series. We more- 
over showed that there was no physical reason to use 
other notions of cosmological distances than the luminos- 
ity distance for evalua ting bounds to the CS, as instead 
previously reported in [29[ . The reason for this lies in the 
fact that dh is adapted the best to the cosmological data 
under consideration. Then, we also showed that the most 
successful parametrization, apart from z, is represented 
by 2/4, as theoretically predicted in Sec. II. 

We carried out our fits up to the sixth order in the CS, 
introducing a further cosmographic parameter, namely 
mo . In addition, we showed that fitting mo together with 
90)ioj s o and lo, is quite feasible to improve the accuracy 
of the analysis by fixing more stringent limits on the CS, 
as opposed to the expectation that an additional fitting 
parameter would significantly broaden the posterior dis- 
tributions. 

The analysis was done for three different sets of pa- 
rameters (including parameters of the CS with different 
maximum order) and two different sets of observational 
data. We essentially found that our numerical results 
appear to be fairly well in agreement with a constant 
pressure associated to the fluid driving the cosmic ac- 
celeration. At first sight, this fluid could be obviously 
thought to originate from a cosmological constant, as de- 
picted in the standard concordance model. Meanwhile, 
our constraints would be able, in future developments, to 
discard different classes of models, which do not satisfy 
the numerical bounds. 

In addition to the numerical fits for the CS, we pro- 
posed a way to constrain both the EoS of the Universe as 
a whole and the derivatives of pressure, by fitting them 
directly from the luminosity distance, up to the order 
of the mo coefficient. In other words, we rewrote dh in 
terms of the EoS parameter and the pressure derivatives 



and performed another MCMC analysis. To achieve this, 
it is necessary to depart from model independent cosmog- 
raphy by assuming a specific gravitational theory. We 
choose General Relativity as valid, and make use of the 
Friedmann equations. The corresponding results lead to 
a set of pressure derivatives compatible with zero within 
the lcr error propagation. Moreover, the EoS parame- 
ter u! is compatible with the one predicted by ACDM,i.c. 
uj = —1/(1 + Slm/f^AO -3 ), obtaining uio — —0.73. Even 
through the pressure derivative results seemed to con- 
firm ACDM, they also left open the possibility that the 
viable model does not necessarily feature a cosmological 
constant, but dark energy with constant pressure and a 
varying barotropic factor, i.e. the VSSM model (as pro- 
posed in [H],[3(|). In this model, which does not involve a 
cosmological constant, the predicted bounds are in agree- 
ment with our fitting results. Unfortunately, the strong 
degeneracy with ACDM in fitting data leaves open the 
question of which model is theoretically the favored one. 

Future perspectives in this direction include using 
more accurate datasets and constraining the analysis at 
higher redshifts, in order to obtain limits able to dis- 
criminate among the two paradigms. We expect that 
this could give further insight into the issue of inferring 
the correct cosmological paradigm, employing a model- 
independent procedure, which does not need to specify a 
model a priori. 



Acknowledgments 

A. A. acknowledges CONACYT for grant no. 215819. 
C.G. is supported by the Erasmus Mundus Joint Doc- 
torate Program by Grant Number 2010-1816 from the 
EACEA of the European Commission. O.L. wants to 
thank Prof. A. Melchiorri for useful discussions. 



[1] A. G. Riess et al., Astrophys. J. 116, 1009-1038 (1998); [2] Garnavich, P. M. et al., ApJ, 509, 74 (1998); R. Rebolo 
S. Perlmutter et al, Astrophys. J. 517, 565 (1999). et al, MNRAS 353, 747 (2004); A. C. Pope et al, As- 



13 



trophys. J. 607, 655 (2004). 
[3] M. Tegmark et al. (SDSS Collaboration), Phys. Rev. D 

74, 123507 (2006); W. J. Percival, S. Cole, D. J. Eisen- 

stein, R. C. Nichol, J. A. Peacock, A. C. Pope, and A. S. 

Szalay, MNRAS 381, 1053 (2007). 
[4] S. M. Carroll, AIP Conf. Proc. 743, 16 (2005); R. Aman- 

ullah et al, Astrophys. J. 716, 712 (2010). 
[5] V. Salmi and A. Starobinsky, Int. J. Mod. Phys. D 15, 

2105 (2006). 

[6] C. Clarckson, G. Ellis, J. Larena, and O. Umeh, Rep. 

Prog. Phys. 74, 112901 (2011). 
[7] E. W. Kolb, S. Matarrese, A. Notari, and A. Riotto, 

Phys. Rev. D 71, 023524 (2005). 
[8] S. Rasanen, JCAP 0611, 003 (2006). 
[9] C. Bonvin, R. Durrer, and M. Kunz, Phys. Rev. Lett., 
96, 191302 (2006). 
[10] A. J. Albrecht et al. , |arXiv:astro-ph /060959T] (2006) . 
[11] A. Torres-Rodriguez and C. M. Cress, MNRAS 376, 1831 

(2007) . 

[12] P. S. Corasaniti, T. Giannantonio, and A. Melchiorri, 

Phys. Rev. D 71, 123521 (2005). 
[13] I. Maor and O. Lahav, JCAP 0507, 003 (2005). 
[14] M. Li, X.-D. Li, S. Wang, and Y. Wang, Commun. Theor. 

Phys. 56 525 (2011). 
[15] J. Einast o, AIP Conf. P roc. 1205, 72-81 (2010); J. 

Einasto, larXiv:0901.0632l (2009); J. L. Feng, Ann. 

Rev. Astron. Astrophys. 48, 495 (2010); LA. Yegorova, 

A. Babic, P. Salucci, K. Spekkens, and A. Pizzella, 

larXiv:1110.1925l (2011). 
[16] A. Aviles and J. L. Cervantes-Cota, Phys. Rev. D. 84, 

083515 (2011). 

[17] M. Berry, Principles of Cosmology and Gravitation, 

Cambdridge University Press, New York (1976). 
[18] T. Padmanabhan, Phys. Rep. 188, 285 (1990). 
[19] S. Weinberg, Cosmology, Oxford Univ. Press, Oxford 

(2008) . 

[20] E. Komatsu e tal, Astrophys. J. Supp. 192, 18 (2011). 

[21] S. Tsuiikawa. larXiv:1004.1493l (2010). 

[22] S. Weinberg, Rev. Mod. Phys. 61, 1 (1989). 

[23] A. Balbi, M. Bruni, and C. Quercellini, Phys. Rev. D 76, 
103519 (2007). 

[24] O. Luongo, H. Quevedo. laxXiv: 1104.47581 (2011). 

[25] B. Ratra and P.J.E. Peebles, Phys. Rev. D 37, 3406 
(1988) R. Caldwell, R. Dave, and P. J. Steinhardt, 
Phys. Rev. Lett. 80, 1582 (1998); G. R. Dvali, G. 
Gabadadze, and M. Porrati, Phys. Lett. B 485, 208 
(2000); C. Armendariz-Picon, V.F. Mukhanov, and P.J. 
Steinhardt, Phys. Rev. Lett. 85, 4438 (2000); R. Micha 
and I. Tkachev, Phys. Rev. D 70, 043538 (2004); E. J. 
Copeland, M. Sami, and S. Tsujikawa, Int. J. Mod. Phys. 



D 15, 1753 (2006); T. Sotiriou and V. Faraoni, Rev. Mod. 

Phys. 82, 451-497, (2010). 
[26] M. Visser, Gen. Rel. Grav. 37, 1541 (2005). 
[27] S. Capozziello, R. Lazkoz, and V. Salzano, Phys. Rev. D 

84, 124061 (2011). 
[28] R. R. Caldwell and M. Kamionkowski, JCAP 0409, 009 

(2004); M. Visser, Class. Quant. Grav. 21, 2603 (2004). 
[29] C. Cattoen and M. Visser, Phys.Rev. D 78, 063501 

(2008) . 

[30] O. Luongo and H. Quevedo, Astroph. and Sp. Sci. 338, 

2, 345-349 (2012). 
[31] C. Cattoen and M. Visser, [arXiv:gr-qc/0703"T22| (2007) . 
[32] O. Luongo, Mod. Phys. Lett. A 26, 20, 1459-1466 (2011). 
[33] L. Xu and Y. Wang, Phys. Lett. B 702, 114 (2011). 
[34] S. Capozziello and V. Salzano, Adv. Astron. 2009, 1 

(2009) ; O. Luongo, G. Iannone, and C. Autieri, Europh. 
Lett. 90, 39001 (2010); G. Iannone and O. Luongo, Eu- 
roph. Lett. 94, 49002 (2011). 

[35] A. G Riess et al., Astrophys. J. 699, 539 (2009); R. 

Jimenez, L. Verde, T. Treu, and D. Stern, Astrophys. 

J. 593, 622-629 (2003). 
[36] A. Lewis and S. Bridle, Phys. Rev. D 66, 103511 (2002). 
[37] T. Faber and M. Visser, Mon. Not. Roy. Astron. Soc. 

372, 136 (2006) 
[38] L. Xu and Y. Wang, Phys. Lett. B 702, 114-120 (2011). 
[39] R. Nair, S. Jhingan, and D. Jain, JCAP 1105, 023 

(2011); B. A. Bassett and M. Kunz, Phys. Rev. D 69, 

101305 (2004). 

[40] I. M. H. Etherington, Gen. Rel. Grav. 39, 1055 (2007); R. 

Kristian and R. K. Sachs, Astrophys. J. 143, 379 (1966). 
[41] J. E. Taylor et al, Astrophys. J. 749, 127 (2012). 
[42] R. Amanullah et al. (SCP Collaboration), Astrophys. J. 

716, 712 (2010). 
[43] R. Lazkoz, S. Nesseris, and L. Perivolaropoulos, JCAP 

0807, 01 (2008); A. Avgoustidis, L. Verde, and R. 

Jimenez, JCAP 0906, 012 (2009). 
[44] R. F. L. Holanda, J.A.S. Lima, and M. B. Ribeiro, As- 
tron. Astrophys. 528, L14 (2011); F. De Bernardis, E. 

Giusarma, and A. Melchiorri, Int. J. Mod. Phys. D 15, 

759 (2006). 

[45] G. F. R. Ellis, Gen. Rel. Grav. 39, 1047 (2007). 

[46] A.G. Riess et al., Astrophys. J. 699, 539 (2009). 

[47] D. Stern, R. Jimenez, L. Verde, M. Kamionkowski, and 
S. A. Stanford, JCAP 1002, 008 (2010). 

[48] J.-Q. Xia, V. Vitagliano, S. Liberati, and M. Viel, Phys. 
Rev. D 85, 043520 (2012). S. Capozziello, R. Lazkoz, and 
V. Salzano, Phys. Rev. D 84, 124061 (2011). 

[49] B. E. Scha efer, Astrophys. J . 660, 16 (2007). 

[50] R. Durrer, larXiv: 1103.5331 1 (2011). 



14 



Appendix A: 5 Distances in terms of redshifts 
1. Redshift z 

We begin with the results for the cosmological distances for the conventional redshift z, starting with luminosity 
distance, 



d L = 



1 

+ z 4 



go 
2 



+ z 6 -( 
, 7j g 3 , 



1 

2 

1 5j gp 

12 24 12 
1 9jo 

20 40 ' 12 

1 7j 19j 2 

30 30 72 

133«g 21g 5 



z • -— - 



Jo 2 



4 48 
the photon flux distance, 



16 



5jogp _ 

12 
_Jp_ + 
120 

720 
13sp 
90 



5g 2 



jo 9o ' 
6 6 
5gp sq 
8 8 24 
go _ Hjoffo 
20 12 
too _ go 4 
720 



27g| 
40 
13j <?o 



30 9 
7jo«o . 19g so 



144 



48 



7jZo£o\ 
24 / 



7j go 


H% 3 


, 7?o 4 


8 


8 


8 


7jo 2 g 


7/o<7o 


7<7o 2 


18 


240 


10 



llso q s 



) 



120 8 ) 
133j gg 13gg 
48 6 



dp = tt- 



1 

Ho 



z 



+ z u 



17?o 



2 

l + ^Zo _ 
24 24 48 

71 47j 

1920 120 12 



_L_ J° 
24 6 

5 jo 9o 



5go 
12 



g 2 



12 



i 2 



1 960 



457jo 11 jo 2 



960 



36 



120 

llZo 

360 



593g 3 , 7j g 3 77g 1 



192 4 
the photon count distance, 



24 



21gp 5 
16 



7? 2 


5<7o , s o N 


8 


8 24/ 


t 149g 9j g 


480 


8 




1069g 


720 


3840 


593s 


7joso 


2880 


144 



, 47gg 
40 
593j g 
288 

llgpSQ 

24 



7jogp 27gg 
8 16 



7j'o9o . 7Z g 



18 



7ggso 
24 



7go ys o go«o \ 
'—J 
77jo9o 



240 



9_so 
80 
457g 2 



320 



24 



dp = — — 



1 



3j 



, 7jogo 



1 

1 

3jo 
5 

I + ^£ 
6 6 

175q 4 



1 _ go 

2 2 

_ 

4 

Jo 2 



/I 3jo jj lp 4g 
V5 5 12 120 5 



i _io 

3 6 

_9gf _ 
12 8 8 

4gp 4j g 

5 



2go 
3 



5 jo go 



5g 3 



24 

9go _ 7jpgg 
5 



8 



4 48 
the deceleration distance, 



72 
21g| 
16 



5Zp 
144 
5so 
18 " 



toq _ 
" 720 " 
7jo5 
144 



5go 
6 



25j go 



+ 2g, 
7jo 2 9o 



9 18 
25g s , 7gggo V 



48 



24 I 



7gp _ 2fo _ go£o\ 



8 
240 



2so 
15 

2 



8 / 
175iog§ 25g 3 



48 



6 



1 

H~ 



z + 



+ z*- - — 



11 
12 

,5 . ( 563 



1 



9o- 
2 

11 jo 61 go 



24 
17 j 



48 



'23 
,24 
5 jo go 
12 



llgo 
12 



30 
6 

Hgo 2 5g 3 



V640 20 
1627 



I 
12 



120 

1273j 7jo 2 



253g 37j gp 



9l 



£0 
24 

51go 



7jog 2 



1920 



960 



18 



160 24 20 8 

7/p too_ _ 7141g 1037j o go 
180 



37g 3 
16 
7j'ogo 



720 



3840 



7g 4 



37s go«o 



8 240 h 
7Z go 1273g 2 



49j g 2 1037g3 7j «§ 49g 1 21g5 1037s 



288 18 240 

7 jo s 7g s . 7g^s 



320 



12 



192 



12 



16 



2880 



144 



12 



24 



15 

and the angular diameter distance, 

4 / 25 13 jo 23go 5joq 13gg 5g^ s \ 

' V 12 24 12 12 8 8 24/ 

+ 5 f 137 _ 137 Jo , _ Jo_ 163g _ 7 jo g 137gg _ 7j gg 21gg 7gp 7s g s \ 

V 60 120 12 120 60 4 40 8 8 8 40 " 8 / 

6 ( 49 79 -jo 31 Jo 31? " m o 71 g" 163 ^gQ 7 ^o 2 go 7l o1o 237( ?o 217 ^o 

V 20 40 72 720 720 20 36 18 240 40 48 
_ 163gg 7j gg _ 21 7q$ _ 21q% 163s _ 7j s 31g s 7g§s y 

24 4 48 16 360 144 48 24 ). ' 



2. Redshift yi 

Here we give the results for the distances in terms of yi, first for the luminosity distance, 
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The cosmological distances are given in terms of yi - luminosity distance, 
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4. Redshift j/3 

These are the results for the cosmological distances in terms of the third redshift, j/3, beginning with the luminosity 
distance, 
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5. Redshift j/4 



Finally, for the last redshift 2/4 the cosmological distances are given, starting with the luminosity distance, 
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Appendix B: The Hubble parameter as a function of redshifts 

We start with the parametrization of the Hubble parameter in terms of the commonly used redshift z, given by 



H(z) 
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Jo ~ % + 



-3go -3gg+ j (3 + 4g ) + s ) + 
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+ 157o (4 + 16 9o + 25^ + Uq* - s ) + 60s + l05q Q s Q + 60 9 2 s ) . 
Further we calculated the expression H(yi), as 
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and in terms of the third redshift, -ff (7/4), 
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Appendix C: Pressure and derivatives as function of redshifts 



In this section we give the expressions for the pressure in terms of the three redshifts, and also its derivatives with 
respect to these redshifts, evaluated at present cosmic time to, which is equivalent to z = yl = 2/4 = 0. 
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The result for the pressure in terms of z reads 
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and finally with respect to 2/4, the results for the pressure and its derivatives are 
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Appendix D: Luminosity distance as function of the EoS parametrization 



These are the expressions for the luminosity distance dz, in terms not of the CS, but in terms of the pressure and 

dP 



its derivatives with respect to redshift. To make the equations more readable, we introduce the notions P\ = 4^, 



P-2 = ^ etc., where y l = z, 3/1,2/4- 
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For rcdshift z we have 
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for redshift y\ the luminosity distance is given by 
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and for the last redshift 2/4 the result reads 
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Appendix E: Derivatives of G(z) 



Here we report the results for the derivatives of the function G(z), which characterizes the equation of state of dark 
energy in a specific model, evaluated at present time, for the three redshifts under consideration. 
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The derivatives with respect to rcdshift z at z = read 
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